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Abstract. The two component plasma (TCP) living in a Flamm's paraboloid is 
studied at a value of the coupling constant F = 2 for which an analytic expression 
for the grand canonical partition function is available. Two cases are considered, 
the plasma in the half surface with an insulating horizon and the plasma in 
the whole surface. The Green's function equation necessary to determine the 
n-particle truncated correlation functions is explicitly found. In both cases this 
proves too complicated to be solved analytically. So we present the method of 
solution reducing the problem to finding the two linearly independent solutions of 
a linear homogeneous second order ordinary differential equation with polynomial 
coefficients of high degrees. In the flat limit one recovers the solution for the 
structure of the TCP in a plane in the first case but the collapse of opposite charges 
at the horizon makes the structure of the plasma physically not well defined in 
the second case. 
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1. Introduction 

The two component plasma (TCP) is a neutral mixture of point wise particles of 
charge ±g. The equation of state of the TCP living in a plane is known since the work 
of Salzberg and Prager [1]. In the plasma the attraction between oppositely charged 
particles competes with the thermal motion and makes the partition function of the 
finite system diverge when P = > 2, where (3 = l/Zc^T with ks Boltzmann 
constant. The system becomes unstable against the collapse of pairs of oppositely 
charged particles, and as a consequence all thermodynamic quantities diverge, so 
that the point particle model is well behaved only for P < 2 [2] when the Boltzmann 
factor for unlike particles is integrable at small separations of the charges. In this 
case rescaling the particles coordinates so as to stay in the unit disk one easily proves 
that the grand canonical partition function is a function of ^^/\ZX^V'-^~'^^'^\ where 
V is the volume of the plasma and X± the fugacities of the two charge species, and 
as a consequence the equation of state is PP = n(l — P/4) where n = n+ + n_ is the 
total particle number density. It also follows that the ratio ^yX^X^/n^^^'^^'^'^ must 
be a function of P only in the thermodynamic limit [3] . However, if the collapse is 
avoided by some short range repulsion (hard cores for instance), the model remains 
well defined for lower temperatures. Then, for P > 4 the long range Coulomb 
attraction binds positive and negative particles in pairs of finite polarizability. Thus, 
at some critical value Pc ~ 4 the system undergoes the Kosterlitz-Thouless transition 
[1] between a high temperature (P < 4) conductive phase and a low temperature 
(P > 4) dielectric phase. The same behavior also occurs in the TCP living in one 
dimension [5]. 

The structure of the TCP living in a plane at the special value P = 2 of the 
coupling constant is solvable exactly analytically [6] . Through the use of an external 
potential it has also been studied in various confined geometries [3, [HI El [TD] and in 
a gravitational field [TTl [12] . 

On curved surfaces it has only been studied in surfaces of constant curvature 
as the sphere [131 E] the pseudosphere [15]. Unlike the one component 
plasma where the properties of the Vandermonde determinant allowed the analytical 
solution a Cauchy identity is used for the solution of the TCP. Unlike in the one 
component case where the solution was possible for the plasma confined in a region 
of the surface now this is not possible, anymore, without the use of an external 
potential. In these cases the external potential is rather given by — (P/g^)lny^ 
where g is the determinant of the metric tensor. 

In this work we study the TCP at the special value P = 2 of the coupling 
constant in a Flamm's paraboloid, a surface of non constant curvature. In this 
surface the one component plasma [16] has been studied in various situations: 
confined to one "universe" by an insulating or a grounded perfect conductor 
"horizon" or in the whole surface. Here we will restrict ourselves only to the first 
and last situation. 
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On a curved surface, even though the finite system partition function will still 
be finite for F < 2 since the surface is locally fiat, the structure will change respect 
to the fiat case. The purpose of the present work is to see how it changes for the 
special case F = 2 where an exact analytical solution is possible. 



2. The Flamm's paraboloid S 

In this work, we want to study a two dimensional TCP in the Flamm's paraboloid 
This is a Riemannian surface S with the following metric 



ds' = gf.^dx^dx" = (^1 - ^ J dr^ + r^dip^ ■ (1) 

or Qrr = 1/(1 - 2M/r),g^^ = r^, and Qr^ = 0. 

The Flamm's paraboloid is an embeddable surface in the three-dimensional 
Euclidean space. It is composed by two identical "universes": 5+ and iS_. These 
are both multiply connected surfaces connected by a hole of radius 2M. We will from 
now on call the r = 2M region of the surface its "horizon" . The scalar curvature is 
R = -2M/r^. 

The system of coordinates (r, ip) with the metric ([1]) has the disadvantage that 
it requires two charts to cover the whole surface S. It can be more convenient to 
use the variable 



u = ±J—-l (2) 
V 2M ^ ^ 

instead of r. This gives the following metric when using the system of coordinates 

(m,V3), 

ds^ = 4M2(1 + u^) [4 du^ + (1 + u^) dif^] . (3) 

The region n > corresponds to 5+ and the region u < to 5_. 
The Laplace- Beltrami operator is 

/, (4) 



where q = (r, ^p). Finding the Green's function of the Laplacian, naturally [16] leads 
to consider the system of coordinates (s,</9), with 

s = {yAFTl + uf. (5) 

The range for the variable s is ]0, +oo[. The lower paraboloid iS_ corresponds to the 
region < s < 1 and the upper one 5+ to the region s > 1. A point in the upper 
paraboloid with coordinate (s, ip) has a mirror image by refiection {u — > —u) in the 
lower paraboloid, with coordinates (l/s,^?), since if 



{V^FTi + uf (6) 
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then 

In the upper paraboloid S+, the new coordinate s can be expressed in terms of the 
original one, r, as 



2M • 

Using this system of coordinates, the metric takes the form of a fiat metric 
multiplied by a conformal factor 



rfs2 = l^h + ^) {ds' + sUif') . (9) 



The Laplacian also takes a simple form 

s 



A/ =^ Aflat/ (10) 



where 



^^.,f = dlf + -dsf+-,dlf (11) 



is the Laplacian of the flat Euclidean space M^. The square root of the determinant 
of the metric is now given by ^ = (M/2)^s(l + s~^)^. 

It is useful to keep in mind the following small M behaviors 

v^= — + — + ^ + 0(M^) . (13) 
3. Coulomb potential created by a point charge 



We here summarize the results found in In Ref. [16] on the determination of the 
Coulomb potential between two charges living in half Flamm's paraboloid with an 
insulating horizon and between two charges living in the full Flamm's paraboloid. 

The Coulomb potential G{s, ip; sq, (fio) created at (s, ip) by a unit charge at 
(so, ipo) is given by the Green's function of the Laplacian 

AG(s, V?; So, v^o) = -27r6^^\s, (p; sq, v^o) (14) 
with appropriate boundary conditions. The Dirac distribution is given by 

6^'^\s, ip; So, ipo) = -^S{s - So)6{<p - <po) (15) 
V9 

Notice that using the system of coordinates (s, (p) the Laplacian Green's 
function equation takes the simple form 

AflatG'(s, p>] So, ipo) = -27r- 6{s - so)6{<p - ipo) (16) 

s 

which is formally the same Laplacian Green's function equation for a fiat space. 
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The Laplacian Green's function equation f ll4p can be solved, as usual, by using 
the decomposition as a Fourier series. Since Eq. f lT^ reduces to the flat Laplacian 
Green's function equation ( lT6l) . the solution is the standard one 

G{s,ip;so,(po) = y^- ( — ) cos [n{(p - If o)]+ go{s, So) (17) 

n=l \ >/ 



where s> = max(s, Sq) and s< = min(s, Sq). 

We consider now the case of particles restricted to live in the half surface (hs) 
iS+, s > 1, by a hard wall located at the "horizon", s = 1. The region iS_, s < 1, is 
empty and has the same dielectric constant as the upper region. We want to consider 
a model where the interaction potential reduces to the fiat Coulomb potential in the 
fiat limit M — )■ 0. The solution of the Laplacian Green's function equation is given in 
Fourier series by Eq. (fT7|) . The 0th order Fourier component go can be determined by 
the requirement that, in the limit M — >■ 0, the solution reduces to the fiat Coulomb 
potential 



I r — r" 



G^^*(r,r') = -In^-^ {U ^ 

Ij 

where L is an arbitrary constant length. Recalling that s ~ 2r/M, when M — 0, 
we find 

M 

5-0(5, So) = -lns> - In— (19) 
and in the half surface 

M 

G "(s, if] So, (/9o) = - In I2; - 2:0! - In — , (20) 

where we defined z = se**^ and Zq = se**^" . And lim^v/ _>.o = G^^^ as desired. 

A similar procedure [TB] gives the Coulomb potential between charges living in 
the whole surface (ws) < s < 00, 

G"^(s, ^- so, ^0) = - In - In ^ , (21) 

with Lo another length scale. 



4. The two component plasma model 

The TCP is a neutral mixture of two species of point charges of charge ±g. In 
this work we want to study the plasma on the Flamm's paraboloid. In Ref. [16] 
the Coulomb potential between two charges living in the whole paraboloid, on half 
surface with an insulating horizon, and in half surface with a grounded horizon was 
found. The first and last case have been summarized in the previous section as they 
will be the two cases considered next. 
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5. TCP in the half surface with an insulating horizon 



When the TCP lives in the half surface with an insulating horizon the Coulomb 



potential is given by Eq. f l20|) . We will use Ui = SiC^"-^^ and Vj = Sje^'^^ to denote the 
complex coordinates of the positively and negatively charged particles respectively, 
where s = (i/r + — 2M)^/2M > 1. Then, following Ref. [6], we use a Cauchy 
identity [17] to rewrite e~^'^'^^ , where U2N is the potential energy of a neutral system 
of positive and N negative charg follows 



exp 



m 


In - 






2Ly^ 


n. 


M ) 




M ) 


det 



+ In 



Vi 



2L/M ' 2L/M 

{Ui - U,j){Vi - Vj 



In 



Ui 



^'3 



2L/M 



Ui 



Ui 



(22) 



(»j)e{i,...,7V}2 

where we had to choose the particular case of a coupling constant F = /3g^ = 2. 

Following Ref. [6], to avoid divergences we start from a discretized model in 
which two sub-lattices U and V are introduced . The positive (negative) particles 
sit on the sub-lattice U(y). Each lattice site is occupied by no or one particle. The 
grand canonical partition function defined as a sum including only neutral systems 
is then 



E 

N=l 



2L 

M 



2N 



N\N 



E 



Ul, 



.,ujsieu 



det 



Ui 



3/ {jj)6{l.-,Af}^ 



(23) 



where the sums are defined with the prescription that configurations which differ 
only by a permutation of identical particles are counted only once, and \ is the 
constant fugacity of the charges of species p = ±. This grand partition function is 
the determinant of an anti Hermitian matrix M explicitly shown in Ref. [7]. 

When passing to the continuum limit in the element Mj^ one should replace Ui 
or Vi by z and Uj or Vj by 2;', i.e. i ^ z and j z' . Each lattice site is characterized 
by its complex coordinate z and an isospinor which is J ) if the site belongs to the 
positive sub-lattice U and I,' ) if it belongs to the negative sub-lattice V. We then 
define a matrix by 

(Jx + icfy'^L/M ax — icry2L/M 



{z\M\z') 



z' 



(24) 



2 z-z' 2 z 

where the a's are the 2x2 Pauli matrices operating in the isospinor space, and 
z = X + iy. 

The matrix A4 can be expressed in terms of a simple Dirac operator 

2L 



{z\M\z') 



M 



aydy] 



In \z 



(25) 
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and the grand partition function can be rewritten as 



{z\M\z') 

= det[I + K~^] , (26) 
with I = 15(2) z') and 

A =A,i±^ + A_i^, (27) 

= XM . (28) 

(29) 

Then, since A^^^ In \z\ = 27iS{s)S{(p)/ s = 2'7iS{z) (where s = \z\ and S{z) is the 
flat Dirac deha function), the inverse operator K is 

K = m~\z)0 , (30) 

where 

m(2;) = m+(2)i^^ + m_(2)^-^ , (31) 
2 

O = j^{(^xd^ + (Tydy) . (32) 

Here m±{z) = {2nLX±y/g/sS){2/MY are rescaled position dependent fugacities, 
g = det((7^i,), and 5* is the area per lattice site which appears when the discrete sums 
are replaced by integrals. Notice that in the flat limit M — we find m± — )■ •rh±, 
where m± = 2nLX±/S are the rescaled fugacities of the flat system [3 [6]. Moreover 
is a local property of the surface independent from its curvature. 
We then find 

In S = Tr {In [I + K"^] } , 

and the one-body densities and n-body truncated densities can be obtained in the 
usual way by taking functional derivatives of the logarithm of the grand partition 
function with respect to the fugacities \±. Marking the sign of the particle charge 
at Zi by an index pi = ±1, and defining the matrix 

Rp,pM, Z2) = {zipi\K~\l + K-y\z2P2) , (33) 

it can then be shown [HI [7] that they are given by 

PpiH^i) = Rp,pAzi,zi) , (34) 

Ppii^(^l'^2) = —Rpip2{Zl, Z2)Rp2pi{z2, Zi) , (35) 
PSI,...,p„(^l'2;2,---,2:n) = 

{il,i2,...,in) 

where mp{z) = {2TiL\p^/ sS){2/Mf = {rhp^ / s){2 / Mf , rhp = 27rLXp/S, ^ = 
(M/2)2s(l + 1/s)^, and the summation runs over all cycles («i, • • • , "^n) built with 
{1,2,. 
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5.1. Symmetries of the Green's function R 
Since mUz) = m(z) and = — O we find 



Rp,p,{zuZ2) = {z2P2\m~'0{I-m-'Oy\in-'0)-^\zm) ■ (37) 
Expanding in O (by using (I— m^^O)"^ = I+m^^O+(m^^O)^ + . . .) and comparing 
with the definition Rp-^p^^zi, Z2) = {zipi\{m^^O)^^{I + m^^O)^^m^^O\z2P2) we find 
Rpp{zi, Z2) = Rpp{z2,zi) , (38) 
Rp-p{zi,Z2) = -R-pp{z2,zi) . (39) 

From which also follows that Rpp{zi, zi) has to be real. If A-|- = A_ then we 
additionally must have 

Rpp{zi, Z2) = R-p-p{zi, Z2) . (40) 

5.2. Two-body truncated correlation functions and perfect screening sum rule 

For the two-body truncated correlation functions of Eq. fl35p we then find 

p?f(^i,^2)= -\R++{ZuZ2)\\ (41) 
P^^I''{Z„Z2) = \R+-{Z,,Z2)\' . (42) 

Notice that the total correlation function for the like particles h^^{zi, Z2) = 

('2)T' ('1') (1) 

pY+ (^15-2^2) / p+ {zi)p\_ {Z2) goes to —1 when the particles coincide zi — Z2 as 
follows from the structure of Eqs. ( !Ml) -( !35l) . Moreover the truncated densities of 
any order has to decay to zero as two groups of particles are infinitely separated. In 
particular \R^+{zi, Z2)\ = |-R++(sl, s2; ip)\, with ip = ipi — (p2, has to decay to zero 
as |si — S2I — 00. 

The perfect screening sum rule has to be satisfied for the symmetric mixture 
j [p+^-izi, Z2) - p+r(^i' Z2)]y/gldsidifi = p±{z2) , (43) 
where gi is g calculated on particle 1. 

5.3. From the structure to the thermodynamics 

Following Ref. [12], we note that if we solve the eigenvalue problem K^^^ = A„'0„ 
(here ipn is a spinor of components tp+^n and ip-^n), then 

^-M'^i)^ (44) 

Now the eigenvalue problem reads 

{m+(^) + A{m^{z)\-^A^}^.[z) = (m„(z) - X'^)tp4z) , (45) 

{m_{z) + A^m+{z)]^^A} ip+{z) = {m_{z) - X'^)ip+{z) , (46) 

where A = {2/M){dx + idy) and A^" = {2/M){-d^ + idy) and the left hand side is 
the same as the left hand side for the Green's functions Eqs. fH^ - (lSU]) . 

Of course we know that the thermodynamic quantities all diverge in our F = 2 
case so Eq. (IHI) . although suggestive, is of small practical interest. 
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5.4- Determination of the Green's function R 

The Green's function matrix R is the solution of a system of four coupled partial 
differential equations, namely 

(I + K-i)KR(zi,Z2) = (I + K)R(zi,Z2) = lS^'\z,;z2) (47) 

where ^(^^(zi; Z2) = 5^'^\s,(f; so,ipo) = {y^)~^5{s — so)5{(p — y^o) is the Dirac delta 
function on the curved surface, 5(s — So)5(v5 — </5o)/s = v?; Sq, v^o) = (^(^S ^0) is 
the fiat Dirac delta function, and 1 is the identity matrix. These can be rewritten 
as follows 

[0 + m(zi)]R(zi,^2) = m(^i)5(2)(^^;^2) . 
If instead of R one uses R = Gm, G satisfies the equation 

[O + m(zi)]G(zi, Z2) = ^2) , (48) 



where S{^Z\ \ Z2) is now the fiat Dirac delta function 

By combining the components of this equation one obtains decoupled equations 
for 6*++ and G as follows 

{m+{zi) + A^[m4zi)]~^A} G++{zi, Z2) = ^^(^i, ^2, ^2) , (49) 
{m_{zi) + A[m+{zi)]~^A^}G^.{zi,Z2) = ^6{si,^i; 82,^2) , (50) 



while 



G^+{zi,Z2)= -[m_(2i)] ^ AG++{zi,Z2) , (51) 
G+_(zi,Z2) = + [m+(^i)]-MtG„„(zi,Z2) , (52) 
Then Eq. (H9|) can be rewritten in Cartesian coordinates as 

4(-xi + iyi) 



siil + s,) i9..+^^.)\\G,AzuZ2] 
2 Y 



V^i + yl 

where s = -^/a;^ + y"^. From the expression of the gradient in polar coordinates 
follows 

dx = COS Lfds , 

dy = sm ipds H d^ . 
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+ + 



4 



+ Si 

A'' 
M 



-d 



Ai 



-9. 



^1 



(55) 



From this equation we immediately see that G-^+{zi, Z2) cannot be real. Notice that 
in the fiat limit M — )■ we have s ~ 2r/M and Eq. fl55l) reduces to 

1 . , . . 1 



—dr,{ridr,) - —d^^ 



ri 



m 



5{ri - r2)S{ifi - (P2) . 



which, when m. 



m_ 



G++{zi, Z2) 



m, has the following well known solution [71 16] 

fh 



27r 



Ko{m\ri - rsl 



(56) 



(57) 



where Ko is a modified Bessel function. 

Let us from now on restrict to the case of equal fugacities of the two species. 
Then A_ = A+ = A with 

^ 27iL , 2^lLe^^' 
m = — —A = 



,2M/g2 



(58) 



S " A2 

where h is Planck's constant, m is the mass of the particles, and fi the chemical 
potential. So m has the dimensions of an inverse length. From the symmetry of 
the problem we can say that (7++ = ^++(51, S2] — ^2)- We can then express the 
Green's function as the following Fourier series expansion 



G 



^ 00 

i,S2]^) = — g++{si,S2]k)e'^^ 
Zn ^ — ' 



(59) 



fc=— 00 



Then, using the expansion of the Dirac delta function, X^fc ^^'^'^ ~ 27r(5((/?), we find 
that g++, a continuous real function symmetric under exchange of si and S2, has to 
satisfy the following equation 

[Qo (/c, si) + Qi (si) ds^ + Q2 (si) 9,^J g++ (si, ss; k) = 



] ms\{l + s{)H (si - S2) 



(60) 



where 



Qo{k,s) =m^{l + s) 



M 



A;s''(4 + A;(l + s)) , 



Qi {s) 

Q2{S) 



2 

M 



/ (5 
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And the coefficients Qi are polynomials of up to degree 9. 



5. 5. Method of solution 

We start from the homogeneous form of Eq. ( 1601) . We note that, for a given fc, the 
two linearly independent solutions /q,(s; k) and //^(s; k) of this linear homogeneous 
second order ordinary differential equation are not available in the mathematical 
literature to the best of our knowledge. Assuming we knew those solutions we 
would then find the Green's function, g+^{si, S2; k), writing [19] 

/(ti, ta; k) = Ckfais<; k)fp{s>] k) , (61) 

where s< = min(si,S2), s> = max(si,S2), and has the correct behavior at large 
s. Then we determine by imposing the kink in / due to the Dirac delta function 
at Si = S2 as follows 

dsJisi, S2; k)\s^=s2+e - dsj{si, S2; k)\s,=s2-e = -m^^^^^^ , (62) 

where e is small and positive. 

The Green's function, symmetric under exchange of Si and S2, is reconstructed 
as follows 



00 



G++{z^,Z2) = G++isuS2■,^) = ^ V cMs^- k)ff,{sy, k)e''^ (63) 

fc=— 00 

6. TCP in the v^hole surface 

On the whole surface, using Eq. fl2T]) . we can now write e~^^^^ at a coupling constant 
r = 2 as follows 



det 



\UjVj\ 



Lq Ui - Vj 



(64) 



(i,i)e{i,...,7V}2 
The grand partition function will then be 

S = det [I + Kws"'] , (65) 

with 

Kws^' = (^] \wsM , (66) 



Ws = — \- X-\z\ — - — . (67) 



The inverse operator is then given by 



Kws = j M-'X^', , (68) 

X+\z\ 2 A_|z| 2 ' ^ ^ 

1 Ss [mV ^ , ^ 
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Introducing position dependent fugacities [2D] 

27r(L/Lo)Ap^ _ ^ 

— = m.^^— — 



Ss 



'P 1 
s 



where now rrip/Lo — nip, we can rewrite 

_ a^ + tay - 
-i^ws — 2 ~' 

with the operators 



z 



m_{z)\z\ 
z 



+ 



2 
1 



m^[z)\z 
1 



-{d^ -idy) , 



{d^ + idy) . 



m-^-{z)\z\'^ m^{z)\z^ 
Then the equation for the Green's functions are 

(1 - a„a+)i2++(2;i, Z2) = S^'^\zi; Z2) , 
(1 - a+a_)R__{zi, Z2) = 5^'^\zi; Z2) , 

-R-i = — a^R , 

R ^ = — . 

The equation for -R++ in the symmetric mixture case 

2 , 25,, dl -td^,+d[ 
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(71) 

(72) 

(73) 
(74) 

(75) 
(76) 
(77) 
(78) 



IS 



m (21) - — 
^1 



„3 
■'I 



s2 
■'I 



TTl^ ( Z\ ) 

^5(Si - S2)S{(fi - ip2) 



■'I 



5(Si - S2)b[}P\ - ^2) 



From this equation we see that R+j^{zi, Z2) will now be real. 

By expanding Eq. ( 1791) in a Fourier series in the azimuthal angle we now 

[Qo {k, si) + Qi{si)ds, + Q2 (si) a,^J g++ (si, S2; /c) = 

^ _ _ V o 



(79) 



find 



where 



(^y^ m.?(l + Si)'5(Sl-.2) 



^0) 



Qi{s) =2s\ 
Q2 (s) = -s^ . 

And the coefficients Qi are now polynomials of up to degree 8. 

In the flat limit we find, for 6*++ = R++ / ni, the following equation 



2d. 



ni - — 



Vi 



-id^^ + (9, 



V>1 



G++{zi, Z2) 



^5(ri - r2)S{ipi - ip2) 
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We then see that we now do not recover the TCP in the plane [71 E]- This has to 
be expected because in the fiat limit the Flamm's paraboloid reduces to two planes 
connected by the origin. 

After the Fourier expansion of Eq. (!59|) we now get 

[Poik,n) + P^{n)dr, + P2ir,)d^.j9+Arur2;k) = m6in - r^) , (82) 



where 



Po{k,r) 


— j^2^ 


Pi{r) 


2 

1 

r 


P2{r) 



k'^-k-2 



The homogeneous form of this equation admits the following two linearly 
independent solutions 



/2(r;-] 
Mr; 2) 

/2(r;2) 

fi{t;k) 
/2(t; k) 



1/2 1 



2mr) 



k 



[D_i/2((-2)^/4v^r) + 



/^-i/2((-2)i/Vmr)]/2 
D_i/2(z(-2)i/4v/^r)]/2 

VrIv 7~4k+4k-^ /4irnry2) 



> k 



else 



where D^{x) are parabolic cylinder functions and /^(x) are the modified Bessel 
functions of the first kind which diverge as /\/27rx for large x ^ — 1/4|. 

Again we write g+^{ri,r2;k) = Ckfa{j<]k)fp{ry^k) and impose the kink 
condition, 

dr^g++iTi,r2]k)\r^=r2+e - dr^g++{ri,r2]k)\r^=r2-e = -fhr2 , (83) 

to find the Cfc. The Green's function is then reconstructed using Eq. ( l63i) . But we 
immediately see that curiously diverges. Even the structure of the plasma is 

not well defined in this situation. The collapse of opposite charges at the horizon 
shrinking to the origin makes the structure of the plasma physically meaningless. 



7. Conclusions 



An attempt was made to extend the solution of the one component plasma, at the 
special value of the coupling constant F = 2, living in the Flamm's paraboloid [16] to 
the two component one, the TCP. The Flamm's paraboloid is a Riemannian surface 
with non-constant curvature which is asymptotically fiat. Its curvature depends on 
a parameter M (2M being the radius of the horizon) in such way that for M -> 
the surface becomes fiat: two planes connected by the origin. 
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The work of F. Cornu and B. Jancovici [TJ |H] on the TCP in a plane showed 
that at such particular value of the coupling constant, using a Cauchy's identity, 
it is possible to find an analytic solution to the structure of the plasma starting 
from a discretized model to avoid divergences and taking the continuous limit in the 
end. All the n-particle correlation functions with n > 1 are well behaved except the 
particle density which diverges as all the other thermodynamic quantities due to the 
competition between the attraction of oppositely charged particles and the thermal 
motion. This can be prevented by adding an hard core to the particles in order to 
prevent the collapse. 

In this work we repeated the same calculation for the plasma living in the 
Flamm's paraboloid, in the half surface case, s > 1, with an insulating horizon (at 
s = 1) and in the full surface case, s > 0. The solution of the equation determining 
the structure, the Green's function equation, is reduced to the mathematical problem 
of finding the two linearly independent solutions of a linear homogeneous second 
order ordinary differential equation with polynomial coefficients of high degree (9th 
in the half surface case and 8th in the full surface case). To the best of our knowledge, 
an analytic solution of the equations turns out to be out of mathematical reach. 
Nonetheless the many-body problem has been reduced to finding the solutions of a 
simple differential equation. The further step of finding explicitly the alanytical form 
of the solutions from the given differential equations is a matter of mathematical 
syntax and we think does not add much to the present work. 

We discussed the symmetries of the Green's function, the perfect screening 
sum rule, the relationship between the structure and the thermodynamics, and the 
method of solution of the Green's function equation. 

We found that for the plasma living in the half Flamm's paraboloid with an 
insulating horizon the flat limit M — )■ reduces the system to the TCP on the 
plane [TJ [6]. For the plasma living in the full Flamm's paraboloid taking the flat 
limit one does not recover the plasma in the plane and this has to be expected 
since the paraboloid in this limit reduces to two planes connected by a hole at the 
origin. Instead the resulting solution for the structure turns out to be physically 
meaningless. The collapse of opposite charges at the horizon as it shrinks to the 
origin spoils the structure of the plasma. 

We leave as an open problem the one of flnding approximations to the two 
linearly independent solutions of the homogeneous counterpart of Eq. (l60l) . for the 
half surface case, and of Eq. flHOj) . for the whole surface CclSG, clS db function of 
the integer k, the further determination of the Green function in these two cases, 
through Eq. (!63l) . the eventual analysis of the solution and discussion about the 
necessity of the addition of a hard core to the particles (hard "disks" of diameter 
D with 6s < 2D{1 + l/s)~^/M) in order to cure the divergence of the densities pp 
with the consequent determination of the equation of state. These flnal steps are 
necessary in order to have the problem of the structure of the fluid reduced to a 
closed form. In the present work we limited ourselves to what can be said in an 
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analytical form. 
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